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Abstract 

In the framework of the recently proposed QCD based parton model for the cumula- 
tive phenomena in the interactions with nuclei two mechanisms for particle production, 
direct and spectator ones, are analysed. It is shown that due to final state interactions 
the leading terms of the direct mechanism contribution are cancelled and the spectator 
mechanism is the dominant one. It leads to a smaller slope of the cumulative par- 
ticle production rates compared to the slope of the nuclear structure function in the 
cumulative region x > 1, in agreement with the recent experimental data. 



1 Introduction 



We have recently proposed a QCD based parton model for the cumulative phenomena in 
the interactions with nuclei [1]. In the model quark diagrams are summed in the vicinity 
of all intermediate thresholds at which some quarks of the nucleus ("donors") transfer all 
their longitudinal momenta to the distinguished active quark and become soft. Particle 
production proceeds in the model via two different mechanisms, the direct and spectator 
ones, schematically shown in Fig. 1 a,b., respectively. The simpler direct mechanism leads 
to the same x dependence of the production rate Id as for the structure function F2{x) in 
the region a; > 1. It is roughly an exponential in x 

Id{x) ~ F2{x) ~ exp(-6oa;) (1) 

where the slope 6o is determined by the QCD coupling constant and quark mass. The 
spectator mechanism besides involves interactions between partons of the projectile and 
target nucleus. It was shown in [1] that each donor quark has to interact with the projectile. 
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As a result the spectator contribution Ig also behaves in x as an exponential but with a 
different slope 

Is{x) ~ ex.p{-bsx) (2) 

the slope bs depending not only on the QCD coupling constant and quark mass but also on 
the partonic amplitude. With a particular parametrization for the latter chosen in [1] the 
slope bs was found to be somewhat larger than bo, so that the spectator contribution resulted 
smaller than the direct one. Then the total particle production rate behaves in x exactly as 
the structure function. At the time when the paper [1] was written there existed no reliable 
data on the cumulative structure functions except the old data on the deuterium at low 
energies and not really in the deep inelastic region [2] . These deuterium data had the slope 
in X much steeper than particle production data from [3,4]. Assuming the direct mechanism 
for particle production this resulted in very different values for the effective coupling constant 
found in [1] for the two cumulative phenomena. In view of a preliminary character of the 
deuterium data we preferred not to pursue this point any further in [1]. 

However recently a new set of data on the structure function of ^^C in the region 0.8 < 
a; < 1.3 was published [5], which essentially agrees with the x dependence observed on the 
deuterium in [2]. Thus it seems to be confirmed that the structure function and particle 
production have different x dependence in the cumulative region. The slope 6o as observed 
in [2] and [5] is roughly 16, whereas the slope bp of the particle production rate is of the 
order 6^8 [3,4]. 

This circumstance gives us a motivation to reconsider our study of the cumulative phe- 
nomena in the light of the experimental evidence for two different exponentials for the struc- 
ture function and particle production. In the framework of our model it means that it is the 
spectator mechanism that gives the bulk of particle production. 

In this note we first demonstrate that final state interactions, not taken into account in 
[1], cancel the leading terms in the direct contribution, so that it becomes much less than the 
spectator one. Then particle production in the cumulative region indeed goes predominantly 
via the spectator mechanism. To be able to explain the experimental slope we study a 
wider class of parametrizations for the partonic amplitude than in [1]. The spectator slope 
bg results very sensitive to the magnitude of the hadronic diffractive cross-sections. The 
parametrization used in [1], which lead to a rather large value for bg ~ 15, corresponded to 
a zero diffractive cross-section. Raising the ratio of the diffractive to clastic cross-section 
up to 1.5 results in lowering bs down to 7^9, not very far from the experimentally observed 
value 6-j-8. However it does not seem possible to push bs still further down without entering 
into a serious conflict with the data on hadronic cross-sections. 

Note that a more phenomenological attempt to explain smaller slopes for the cumulative 
particle production is made in [6], where the existence of multiquark clusters in the nucleus 
and their properties are postulated and the quark-gluon string model formalism of [7] is used 
to calculate production rates. 

The paper is organized as follows. In Sec. 2 we show how the leading terms in the 
direct contribution are cancelled by final state interactions. In Sec. 3 the partonic amplitude 
is related to the hadronic one and the data on hadronic cross-sections are used to study 
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its possible parametrizations. In Sec. 4 we calculate the particle production rate by the 
spectator mechanism and compare its x dependence to that of the structure function and 
the experimentally observed one. 



2 Cancellation of the direct contribution 

In the upper blob of Fig. 1 a the inclusive cross-section off the active quark enters in the 
triple Regge kinematical region, since the scaling variable x of the produced particle should 
be as close as possible to that of the active quark. Therefore we can redraw the diagram 
of Fig. 1 a in the form shown in Fig. 2 a. However we have also to take into account final 
state interactions of the produced particle with donor quarks which raise its longitudinal 
momentum. An example of these is shown in Fig. 2 b. The diagrams of Fig. 2 a and b 
are not the only ones which give contribution. In fact taking a diagram with a certain 
pattern of interactions with donor quarks one should sum over all places where the two 
upper reggeons from the triple Regge interaction are attached to the active quark. We are 
going to demonstrate that this sum vanishes in the limit when the donor quarks loose all 
their longitudinal momentum. 

Let the number of donors be n. We denote the momenta of the left active quark kq, ...ki 
before the interaction with the reggeon and kJ, k,'i_^i...k,'^ = k, after this interaction. Analogous 
quantities on the right side hereafter will be labelled with tildas (see Fig. 2 a, 6). We use 
the hght-cone variables k — {k±,k±). The donor quarks he on the mass shell between the 
interactions. So there minus momenta ("energies") take on physical values 

ki- = (m^ - kij_)/2ki+ = fii (3) 

From energy conservation we find for the energies of the active quark 



N-l 



and 



Ki_ = Np^_-J2l^j-Ef'j-12kj-^ z = l,2,...Z (4) 

n+l 1 i+1 



= k_ + ^{i^j - kfl) i = 1,1 + 1,. ..n-l (5) 



where N is the total number of partons in the nucleus and kf'^ are the momenta of the 
left donors before the interaction. Integrations over k\^ can then be done easily, since all 
singularities coming from the active quark propagators lie on the opposite side of the real 
axis as compared to the donor propagators. As a result kf^ arc substituted by nf'^ in (4) 
and (5). In the limit when the donor quarks loose all there longitudinal momenta //j become 
large and from (4) and (5) one obtains for the ith active quark propagator 

^=1,2,.../ (6) 
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and 

n 

-1/^/.,-, z = /,/ + l,...n-l (7) 

i+l 

Now we pass to the triple Regge interaction. The momentum of the upper left reggeon 

q is 

q^K'i-Ki (8) 
So Q- — jii. Calculating g+ and q± we find 

= -2pi+(A - a;,) E //^ + + E ^Ol (9) 
11 1 

Here Xi are the scaling momenta of the donors after the interactions, Xi ^ 0; A = n + 1 — x 
where x is the scaling variable of the produced particle. All scaling variables are defined here 
relative to the average longitudinal momentum of the initial quarks of the nucleus. They are 
three times larger than the standard ones defined relative to the nucleon momentum. At the 
threshold A — > 0. The important thing about (9) is that q'^ does not depend on I, that is, 
on the place at which the reggeon is attached. 

However there still remains some dependence on / in the energy on which the upper 
reggeons depend. In fact the factor corresponding to the triple Regge interaction is given by 

t^Q^ (10) 

Here q"^ is given by Eq. (9), g{t) is the three-reggeon vertex, = (A — J2iXi)s where s 
is the standard energy variable; s; = ^is where is the scaling variable of the left active 
quark at the moment of its interaction with the reggeon and sj is a similar quantity on the 
right. All the three reggeons are taken to be pomerons with the trajectory a(t). We have 
also to take into account that the Ith propagator of the active quark is splitted into two by 
the interaction with the reggeon, which results in an extra factor 1/^;. Introducing it into 

(10) we finally obtain for the triple Regge interaction factor 

g{t)s'^^'\^ii!r^'^-'{A -±x,r^'^-"^^'^ (11) 

1 

If we assume weak (logarithmic) dependence of the cross-section on energy then we have 
to take for the effective pomeron intercept a{0) = 1. Then the dependence of the expression 

(11) on / and / is also weak (it enters only through the term * with small a' and t). 
Neglecting this dependence we find that the triple Regge interaction factor does not depend 
on the place where it is attached and is common to all the diagrams of the type shown in 
Fig. 2. 

Then we have only to sum various diagrams forgetting about the common triple Regge 
factor. Diagrams with different / (or/and /) differ only by their active quark propagators. 
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For a given / the factor coming from the left active quark propagators is found from (7) and 

(8) to be 

I i Ji— 1 n 

IKVE/^.OIK-VE/^.) (12) 

1=1 1 i=l i+1 

Upon symmetrizing in the soft donors it becomes 

n 

(n(n-0!)-^(-l)"-^niM (13) 

1 

The sum of (13) over all / = 0, 1, ...n evidently gives zero. Thus in the described approxi- 
mation when the slow change of the triple reggeon interaction with the energy of the active 
quark is neglected all the diagrams of the type shown in Fig. 2 cancel in the sum. 

3 Parton interaction and hadronic cross-sections 

With the bulk of the direct contribution cancelled, particle production in the cumulative 
region goes predominantly via the spectator mechanism of Fig. 1 b. It involves interactions 
between partons of the projectile and donors from the nucleus, characterized by the partonic 
amplitude a. We normalize it according to 

2Ima(0) = (7 (14) 

where a(0) is the forward amplitude and a is the partonic cross-section. As shown in [1] each 
donor parton has to take part in the interaction with the projectile. With n donors, this 
gives a factor a" which is responsible for the difference between the slopes of the structure 
function and particle production. The parametrization of a and its magnitude thus aquire 
the decisive role in describing the experimental slope. 

Of course, a is not the quantity to be directly measured experimentally. It can however 
be related to the data through hadronic {pp or pp) interaction [8]. The elastic hadronic 
amplitude can be represented through partonic interaction as shown in the diagram of Fig. 3. 
Taking the cm. system one finds that the longitudinal components gj± of the transferred 
momenta arc small. Also the upper part of the diagram belonging to the projectile does not 
depend on gj+ and the lower, target part does not depend on Integrating over qi± one 
then finds that for a diagram with n interactions, M partons in the projectile and N ones 
in the target the amplitude is given by 

.Af '^)(g^) = n\Cl,C% I ll{j^,ta{q.±)){27rmq - Fi^\^^±) (15) 

Here F^'^'-''^\qi±) are the n-fold transverse form-factors for the projectile (M) and target 
(N). Their Fourier transforms give parton distributions in the transverse space F^^'^\bi). 
The amplitude (15) thus can be presented as an integral in the transverse space 

/n 
d'BeM^QB) l[{d\d'b'AB -b, + b',))Fi^\h)F^''\b',) (16) 
1 
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This is only a contribution from given numbers of partons in the projectile and target. 
Summing over M and N we obtain for the elastic amplitude with n partonic interactions 

/n 
d'B eMiQB) m^b.d'bXB - 6^ + h'^))F,,{hi)FM (17) 
1 

where for the projectile 

F„(6,)= E(^!/(M-n)!)FW(6,) (18) 

M>n 

and similarly for the target. 

To move further one has to make some assumptions about the multiparton distributions 
Fn{bi). The simplest one is to assume that partons are completely independent: their number 
is distributed according to Poisson's law and multiparton distributions factorize. One then 
obtains 

n 

F^{h,...b^)^u^l[p{h) (19) 
1 

where v is the mean number of partons in the projectile or target and the one parton 
distribution p is normalized to unity 

J d%p{b) = 1 (20) 

Putting (19) into (17) and summing over n one arrives at an cikonal amplitude 

iA{q±) = j d^Bex.p{iqB){exp{ip{B)) - 1) (21) 

where the eikonal factor is 

p(B) = UpUt J d%d%'pp{b)pt{b')a{B -b + b') (22) 

and the subscripts p and t refer to the projectile and target, respectively. This form of the 
amplitude coincides with the one found in the multipomeron exchange model with factoriz- 
able vertices [9]. Eq. (22) then gives the contribution of a single pomeron exchange. 

One commonly assumes that both p and a have a Gaussian dependence on the impact 
parameter: 

= (lMp,t)exp(-rVrp,t), a(r) = (icT/27rr^) exp(-rVro) (23) 

where we have taken a pure imaginary, for simplicity. Then one easily finds simple expressions 
for the eikonal factor, the hadronic total, elastic and inelastic cross-sections and also for the 
slope of the elastic cross-section B^'': 

p{B) = ixexp(-SVi?^) 
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(7'^' = rf(20i(x)-0i(2x)) 
(7^" = 7ri?Vi(2x) 

S^' = (it:V2)02(x)M(x) (24) 

where i?^ = + + Tq is the total interaction radius squared, x — Vpftcr /{2t:B?) and the 
functions 01,2(2^) are defined by 

00 

ct>n{x) = Y.{-lf-'x^/{k\k^) (25) 

k=l 

Taking the ratio of any pair of the quantities in (24) and comparing it to its experimental 
value one can determine x. Then one of Eqs. (24) can be used to find i?^. With x and B? 
determined, the cross-section VpVtcy can be found. Values of Vp and rt are more or less known 
from electromagnetic properties of the projectile and target, thus tq can also be found. So, 
in the end, the only parameters left are the average numbers of partons in the projectile and 
target Vp and Vf At low energies the valence quark approximation seems to be good enough, 
which implies Vp = Vt = ?> for pp or pp interactions. 

This standard procedure was used in [1] to determine the parameters entering the spec- 
tator mechanism. As mentioned, it lead to a relatively small a, so that the slope of the 
spectator spectrum resulted even steeper than that of the structure function. However, 
apart from this unsatisfactory result for cumulative production, the discussed parametriza- 
tion has a more fundamental (and well-known) defect. As one can deduce from the first 
three equations in (24) it gives no diffraction. This property of the pure eikonal amplitude 
can be easily understood in the framework of the Gribov approach to multiple scattering 
[10], where it corresponds to retaining only the initial particle state in the sum over all in- 
termediate states. A possible remedy consists in changing the eikonal amplitude (21) by the 
quasi-eikonal one with a diffraction factor ^ > 1 [11]: 

iA{qi_) = j (fBexi){iqB)^-\exp{i^p{B)) - 1) (26) 

With this factor Eqs. (24) change as follows 

(T*"* = r'2rf0i(x) 

cr"' = r^7rR^{2(f)i{x) - 0i(2x)) 

(7^" = C^TTR^<f)i{2x) (27) 
and S^' does not change at all. Prom (27) one finds the diffractive cross-section 

^d^f _ _ ^yeZ (28) 

SO that the new parameter ^ can be directly taken from experiment. 
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In our partonic approach the quasi-eikonal parametrization (26) means that instead of 
(19) we take 

n 

i^n(&l,-M=e^""'^/'i^"nP(^i) (29) 

1 

both for the target and projectile. 

The quasi-eikonal parametrization leads to values of a considerably larger than without 
diffraction because of a stronger screening effect introduced by diffractive states (sec Table). 
For the maximum value of ^ ~ 2.4 compatible with the experimental data at \/s = 23.5 GeV 
it leads to the value of the corresponding parton-nucleon cross-section ai nearly three times 
larger than with ^ = 1 used in [1]. Note that the absolute value of ai then results quite 
large (~ 50 mb). It is interesting that values of ai of a similar order are also favoured by the 
study of the behaviour of the nuclear structure functions in the opposite kinematical region 
at small x [12]. In relation to this we recall that the parton-nucleon cross-section ai is not 
a directly observable quantity. It only enters hadronic cross-sections as a parameter. The 
mentioned large value of ai is consistent with all experimental evidence on proton-proton 
interactions at ^/s = 23.5 GeV. This large value transforms into the physical value of the 
pp cross-section as a result of a strong screening effect, that is, because of a highly coherent 
manner in which the quarks interact. 

4 Numerical results and discussion 

With the partonic amplitude fixed in the preceding section, one can calculate the spectator 
contribution to the cumulative production rate, using the coresponding equations from our 
paper [1] (Eqs. (38), (39) and (41) of [1]), suitably modified to the quasi-eikonal case. We 
have also corrected a mistake eliminating p\ in the denominator of Eq. (39) of [1]. We have 
chosen the maximum possible value for the diffractive parameter ^ = 2.4. Other parameters 
have been fixed as in [1]. 

The results of the calculations are shown in Figs. 4 and 5 for the cumulative charged pion 
production on deuterium and ^^^Ta, respectively. For comparison the corresponding nuclear 
structure function ^2(2;) at a; > 1, calculared in [1], is also shown, as well as the available 
experimental data from [2-5]. 

One clearly observes that the spectaor mechanism, with a parametrization of the partonic 
amplitude chosen to account for diffraction, leads to a considerably smaller slope of the 
production spectra {bg ~ 7-^9, see Introduction) compared to the slope of the structure 
function in the region a: > 1 (69 ~ 16), in a good agreement with experimental data. 

Thus our model correctly predicts two differents exponential in the cumulative produc- 
tion and in the nuclear structure function at a; > 1. The difference is due to additional 
multiple interactions between projectile and target which enter the spectator mechanism for 
the cumulative production. 
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Table. Cross-sections and elastic slopes for various 
values of the diffract ive factor ^ 







(Ti {mb) 


B^' {GeV-'') 


1.0 


0. 


16.5 


11.1 


1.5 


3.4 


20.2 


10.9 


2.0 


6.8 


31.2 


10.6 


2.2 


8.2 


44.3 


10.6 


2.4 


9.5 


48.3 


11.4 



Table captions 

The first column gives values of the diffractive parameter ^. In the second column the 
corresponding pp diffractive cross-section at ^/s — 23.5 GeV is shown. In the third column 
the values of the parton-projectile (proton) cross section are shown, as calculated from 
Eqs. (24) (from the value of x with z/ = 3). The last column shows the resulting elastic 
slope, its experimental value at ^/s = 23.5 GeV being B^^') = 11.8 ± 0.30 GeV'^ [13]. 
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Figure captions 

Fig. 1 Two mechanisms of the cumulative particle production in the parton model, the 
direct (a) and spectator (b) ones. Dashed and chain lines show gluon and pomeron 
exchanges, respectively. 

Fig. 2 The direct mechanism of the cumulative particle production in the triple Regge 
approach without (a) and with (b) final state interactions of the produced particle 
with donor quarks. Shown for the cases {n — I = I = 2) and (n = 2, / = 1, / = 0), 
respectively. Notations as in Fig. 1. 

Fig. 3 The elastic hadronic amplitude in the parton model. Notations as in Fig. 1. 

Fig. 4 Id = ^(mb) is the calculated inclusive cross-section (per nucleon) for cumulative 
charged pion production on deuteron at ^/s — 23.5 GeV (solid curve) and 1800 GeV 
(dashed one). F-f /2 is the calculated structure function of the deuteron (per nucleon) 
for a; > 1 at = 6 GeV'^ (sohd curve), 20 GeV'^ (dashed curve) and 500 GeV'^ 
(short dash curve). * - the experimental data [2] on the deuteron structure function 
at 0.8 < < 6 GeV'^. The experimental errors are much less than the star symbols 
and not discernable on our scale. 

Fig. 5 lA{rnb) is the same as in Fig. 4 but for the production on ^^^Ta. A - the experimental 
data [4] on the cumulative charged pion production on ^^^Ta by 400 GeV incident 
proton beam. -F^/2 is the calculated nuclear structure function for the ^^^Ta at = 
50 GeV"^ (dashed curve) and ^^C at = 100 GeV^ (sohd curve). □ and x - the 
experimental data [5] on the ^^C structure function at = 61 GeV^ and 150 GeV^, 
respectively. 
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Fig. 4 
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